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ABSTRACT 

A Li-rich red giant star (2M19411367+4003382) recently discovered in the direction of NGC 6819 
belongs to the rare subset of Li-rich stars that have not yet evolved to the luminosity bump, an 
evolutionary stage where models predict Li can be replenished. The currently favored model to 
explain Li enhancement in first-ascent red giants like 2M19411367+4003382 requires deep mixing into 
the stellar interior. Testing this model requires a measurement of 12 C/ 13 C, which is possible to obtain 
from APOGEE spectra. However, the Li-rich star also has abnormal asteroseismic properties that 
call into question its membership in the cluster, even though its radial velocity and location on color- 
magnitude diagrams are consistent with membership. To address these puzzles, we have measured a 
wide array of abundances in the Li-rich star and three comparison stars using spectra taken as part of 
the APOGEE survey to determine the degree of stellar mixing, address the question of membership, 
and measure the surface gravity. We confirm that the Li-rich star is a red giant with the same overall 
chemistry as the other cluster giants. However, its logg is significantly lower, consistent with the 
asteroseismology results and suggestive of a very low mass if the star is indeed a cluster member. 
Regardless of the cluster membership, the 12 C/ 13 C and C/N ratios of the Li-rich star are consistent 
with standard first dredge-up, indicating that Li dilution has already occurred, and inconsistent with 
internal Li enrichment scenarios that require deep mixing. 

Subject headings: open clusters and associations: individual (NGC 6819)-stars: abundances-stars: 
chemically-peculiar—stars: late-type 
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1 . INTRODUCTION 

Anthony-Twarog et al. (2013, hereafter AT13) recently 
reported the discovery of a Li-rich red giant (RG) star 
in NGC 6819. This star, 2M19411367+4003382, has 
A(Li)~ 2.3 dex ([Li/H]~1.3) and is unusual even within 
the class of rare Li-rich giants. The star’s current po¬ 
sition on a color-magnitude diagram (CMD) suggests 
that it has evolved beyond the Li dilution phase of first 
dredge-up (FDU), and, indeed, the other NGC 6819 
RGs at similar magnitudes have A(Li)< 0.7 dex (AT13). 
What makes this Li-rich star so unusual is that its CMD 
position is clearly below and blueward of the luminosity 
bump, the evolutionary stage where models have demon¬ 
strated that newly synthesized Li can be circulated to the 
stellar envelope, resulting in a brief stage of Li-richness 
(see, e.g., Charbonnel & Balachandran 2000, Eggleton 
et al. 2008, Denissenkov 2012, hereafter D12). Other re¬ 
cent studies (e.g., Kumar et al. 2011 and Carlberg et al. 
2012) have identified Li-rich stars in the field that are 
too hot to be bump stars. The former study hypothe¬ 
sized that these stars are red clump (RC) stars that re¬ 
plenished Li during the He flash, and the low 12 C/ 13 C of 
many of these stars are consistent with internal Li regen¬ 
eration. The latter study favored planetary engulfment 
because their stars tend to be more rapidly rotating and 
have 12 C/ 13 C near normal FDU values. 

The currently favored model for 2M19411367+4003382 
is the mixing model of D12, which can explain Li-rich gi¬ 
ants below the luminosity bump. In addition to enriching 
the surface with Li, the mixing alters the star’s evolution, 
making it follow a more extended “retrograde” evolution¬ 
ary path down the red giant branch (RGB) compared to 
normal luminosity bump evolution. In the discovery pa¬ 
per, AT13 argued that 2M19411367+4003382’s Li abun¬ 
dance and position on the RGB was most consistent with 
this model, but they also noted that a crucial test of this 
hypothesis is the measurement of 12 C/ 13 C to confirm 
whether the extra mixing has indeed occurred. 

What makes 2M19411367+4003382 such a powerful 
test of Li-regeneration hypotheses is its apparent mem¬ 
bership in an open cluster, which allows a direct compar¬ 
ison of its abundances to other stars of nearly identical 
evolution. However, its membership is not conclusively 
established. Its radial velocity (RV), original proper mo¬ 
tion measurement (Sanders 1972), color, and magnitude 
(e.g, AT13) support membership (see Figure 1), but its 
asteroseismic properties (Stello et al. 2011) and more re¬ 
cently measured proper motion (Platais et al. 2013) sup¬ 
port non-membership. Both of these non-membership 
criteria have caveats. Platais et al. (2013) reported that 
stars with high RV membership probabilities but low 
proper motion membership probabilities tend to have 
larger proper motion errors that stars of similar bright¬ 
ness. This suggests that there may have been confusion 
in the star’s identification on the photographic plates. 
The asteroseismic properties of the Li-rich star suggest a 
lower log g and mass than expected for cluster member¬ 
ship, but they yield a luminosity that places the star near 
the distance of the cluster. (These caveats are discussed 
in more detail in Section 7.2). 

In this paper, we use high-resolution, infrared spectra 
to address questions on the Li-rich star’s cluster member¬ 
ship and evolutionary status. We select three similarly 
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Fig. 1.— Top: 2MASS color-magnitude diagram of NGC 6819 
member stars (black squares), 2M19411367+4003382 (large blue 
circle), the three comparison stars (small orange circles), and the 
star overlapping the optical study of Bragaglia et al. (2001, pur¬ 
ple asterisk). Bottom: Radial velocity histogram of NGC 6819 
members. The x marks the RV of the Li-rich star. 
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evolved comparison RGs in NGC 6819 (Section 2) and 
measure abundances of ten elements in the Li-rich star 
and comparison stars (Section 3). We measure a spectro¬ 
scopic log < 7 , and find that it is consistent with the astero¬ 
seismic logg (Section 4). The APOGEE spectra also re¬ 
veal that the Li-rich star is rotating slightly more rapidly 
than the other RGs (Section 5). We combine our data 
and literature resources to consider possible contamina¬ 
tion from companions or unrelated background sources 
(Section 6 ). Finally, we explore the implications of these 
results under the two possible cases of cluster member¬ 
ship (Section 7) and summarize our findings (Section 8 ). 

2 . SELECTION OF COMPARISON STARS 

The stars in this paper were observed with the Sloan 
Digital Sky Survey (SDSS) telescope (Gunn et al. 2006) 
and the Apache Point Observatory Galactic Evolution 
Experiment (APOGEE) instrument (Allende Prieto et 
al. 2008; Majewski 2012) as part of the open cluster 
calibration sample for the APOGEE survey (Zasowski 
et al. 2013; Meszaros et al. 2013). We analyzed the 
Li-rich star and three similarly evolved RGs in NGC 
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TABLE 1 

Measured Properties of NGC 6819 Red Giants 


Property 

Li-rich 

Comparison #1 

Comparison #2 

Comparison #3 



(Most Similar) 

(Near RC) 

(Near Bump) 

2MASS 

2M19411367+4003382 

2M19404965+4014313 

2M19412222+4016442 

2M19412176+4012111 

KIC 

4937011 

5111940 

5112744 

5112734 

T e ff (photometric) 

4700 

4725 

4700 

4620 

T eff (ASPCAP) 

4670 

4687 

4645 

4632 

log L/Lq 

1.52 

1.59 

1.67 

1.69 

log g a 

2.8 

2.7 

2.6 

2.6 

£ (km s -1 ) 

1.5 

1.5 

1.5 

1.5 

A(Li) (dex) b 

2.3 

$ 0.7 

$ 0.7 

$ 0.7 



APOGEE- derived 



logs 

2.35 

2.53 

2.51 

2.62 

A( 12 C) 

8.48+0.06 

8.34+0.05 

8.39+0.05 

8.38+0.02 

A( 14 N) 

8.11+0.05 

8.15+0.02 

8.13+0.02 

8.16+0.03 

A( le O) 

8.94+0.1 

8.82+0.05 

8.75+0.04 

8.77+0.02 

C/N 

2.3 

1.6 

1.8 

1.7 

12 C/ 13 C 

25+5 

24+5 

25+5 

20+8 

A(Na) 

5.98+0.14 

6.28+0.18 

6.20+0.04 

6.19+0.09 

A(Mg) 

7.44+0.11 

7.48+0.11 

7.45+0.13 

7.47+0.12 

A(A1) 

6.32+0.03 

6.47+0.05 

6.45+0.04 

6.46+0.04 

A(Si) 

7.50+0.06 

7.51+0.06 

7.51+0.08 

7.53+0.09 

A(Ti) 

4.93+0.17 

5.02+0.07 

5.02+0.04 

5.03+0.12 

A(Fe) 

7.45+0.12 

7.52+0.04 

7.48+0.12 

7.48+0.08 

A(Ni) 

6.22+0.07 

6.19+0.06 

6.22+0.08 

6.22+0.05 

usini (km s —1 ) 

8.5+1.1 

6.6+1.6 

4.4+1.5 

5.8+1.8 



ARUKASC 



^max (/iHz) 

29.51 

52.33 

43.69 

40.11 

A v (/iHz) 

4.15 

5.16 

4.39 

4.12 

logs 

2.37 

2.62 

2.54 

2.50 

R/Rq 

9.1 

10.4 

12.0 

12.4 

M/Mq 

0.7 

1.7 

1.8 

1.8 

log L/Lq 

1.56 

1.69 

1.80 

1.80 


a log g from the isochrone and used to derive abundances. 

k Measurement for 2M19411367+4003382 from AT13; their Figure 3 provides estimates for comparison stars. 


6819, using an infrared CMD to select stars most sim¬ 
ilar to the Li-rich star, as illustrated in Figure 1. All of 
the NGC 6819 member stars observed by APOGEE are 
shown in Figure 1 together with a 2.25 Gyr, Z = 0.023 
isochrone (Marigo et al. 2008) shifted by (m — M) 0 = 
11.85 (Basu et al. 2011), and E(J - K s ) = 0.07 (from 
E (B — V) = 0.14 mag, Bragaglia et al. 2001). The 
Li-rich star sits below the RC and blueward of the lu¬ 
minosity bump of the isochrone. The same is true of 
the most similar control star (2M19404965+4014313). 
Another comparison star is near the red edge of the 
RC (2M19412222+4016442), and the remaining compar¬ 
ison star (M19412176+4012111) overlaps the luminosity 
bump within the uncertainties. The 2MASS and KIC 
identifiers of these fours stars are listed in Table 1. Addi¬ 
tionally, 2M19411476+4011008 was selected for analysis 
solely because it has previously reported spectroscopic 
abundances and log g derived from optical spectra (Bra- 
gaglia et al. 2001, hereafter B01). This star, identified as 
978 in B01, is hotter than the Li-ricli star by ~ 200 K 
and sits at the blue edge of the RC. Because this star 
is much more evolved, it is not a good comparison star 
to the Li-ricli star, and we provide results for this star 
separately in Table 2. 

From the 2MASS colors and isochrones, we obtain the 
temperature, gravity and luminosity listed in the top sec¬ 
tion of Table 1. T e s is the average of two (J—K s ) 0 to T e s 
calibrations (Bessell et al. 1998 and Gonzalez Hernandez 
& Bonifacio 2009). The uncertainty in (J — K s )o from 
the 2MASS photometry is typically 0.028 mag, leading 


TABLE 2 

Properties of 2M19411476+4011008 compared to literature 


Property 

this work 

B01 

A (this work — B01) 

T e fi 

4920 

4855 

+65 

log L/Lq 

1.79 



log S a 

2.6 

2.60 

+0.00 

log S b 

2.58 

2.60 

-0.02 

£ (km s” 1 ) 

1.5 

1.26 

+0.24 

A( 12 C) 

8.32+0.08 



A( 14 N) 

8.21+0.03 



A( 16 0) 

8.74+0.16 



C/N 

1.3 



12 C/ 13 C 

> 20 



A(Na) 

6.46+0.07 

6.96 

-0.50 

A(Mg) 

7.56+0.09 

7.65 

-0.09 

A(A1) 

6.56+0.02 

6.44 

+0.12 

A(Si) 

7.60+0.03 

7.82+0.14 

-0.22 

A(Ti) 

5.04+0.08 

5.22+0.16 

-0.18 

A(Fe) 

7.58+0.06 

7.60+0.12 

-0.02 

A(Ni) 

6.28+0.06 

6.37+0.14 

-0.09 

usin i (km s —1 ) 

3.4+1.3 




a log g from the isochrone, 
k log g from ionization balance. 


to a temperature uncertainty of ~ 80 K in each calibra¬ 
tion. The Bessell et al. (1998) calibration uses the de¬ 
rived color indices for ATLAS9 overshoot models (their 
Table 1) and gives slightly cooler temperatures than the 
Gonzalez Hernandez & Bonifacio (2009) calibrations by 
40-60 K. This systematic difference in temperature is 
well within the ler uncertainties of each calibration. The 
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resulting temperatures are also in good agreement (all 
within 55 K) with those derived by the APOGEE Stel¬ 
lar Parameters and Chemical Abundance Pipeline (AS- 
PCAP, Garcia Perez, et ah, in prep.). We can also tie 
this temperature scale to an independent spectroscopic 
analysis. The comparison star overlapping the B01 study 
has T e ff = 4855 K derived from optical spectra compared 
to Taspcap = 4867 K, and our T e g = 4920 K. 

Assuming a stellar mass of 1.7 M 0 from the isochrone, 
we derive the surface gravities in the top section of Ta¬ 
ble 1. These are the logo’s used to derive abundances. 
For all of the models, we adopted a microturbulence £ of 
1.5 km s -1 , following the prescription adopted by ASP- 
CAP (£ = 2.24- 0.3 log g ) 24 for log g = 2.6 dex.‘ 

3. ABUNDANCES FROM APOGEE SPECTRA 

The spectrum for each star has been processed with the 
APOGEE pipeline (Nidever et al. 2015). The observed 
data used here are the continuum-normalized “aspcap- 
Star” spectra, with the exception of the wavelength re¬ 
gion near the Ti II line, which falls near one of the detec¬ 
tor gaps and is masked out in the “aspcapStar” spectra. 
This spectral region comes from the “apStar” spectra, 
which we continuum normalize. All spectra are described 
and are available in the tenth data release (DR10, Ahn 
et al. 2014) of the SDSS-III (Eisenstein et al. 2011) 25 . 

Abundances are measured via spectrum synthesis us¬ 
ing MOOG (Sneden 1973) 26 and the APOGEE linelist 
(Shetrone et. al., in prep.), focusing on small subsets of 
atomic and molecular features identified for each element 
by Smith et al. (2013, hereafter S13). All of the wave¬ 
lengths in this paper refer to air wavelengths. A spec¬ 
trum of a ~20 A region around each feature under con¬ 
sideration was generated using MOOG, smoothed by a 
Gaussian to match the broadening of the observed spec¬ 
tra. Small adjustments were made to the continuum level 
and velocity shift of the observed data to best match the 
models. 

3.1. C, N, O and 12 C/ 13 C 

The abundances of C, N, and O must be found itera¬ 
tively because they are inter-related by molecular equi¬ 
librium conditions. We first measure C abundances from 
CO lines, then fix the C abundances to measure O using 
OH lines. These steps are iterated until convergence is 
reached. The N abundances can then be derived from CN 
lines. The lines used in this analysis are from Table 4 of 
S13, excluding the (4-l)V-R 12 C 16 0, (3-l)P 2 9.5 16 OH, 
and (0-1) R1 68.5 12 C 14 N lines. These abundances and 
the C/N values are listed in the middle section of Ta¬ 
ble 1 for the Li-rich star and comparison stars and in 
Table 2 for 2M19411476+4011008. The uncertainties are 
the standard deviations of the individual measurements, 
and we address the uncertainties propagated from the 
stellar parameter uncertainties in the next section. The 
three comparison stars all have nearly identical abun¬ 
dances, but the Li-rich star appears to have somewhat 
higher O. The low C/N confirms that standard FDU has 
completed and Li dilution should have occurred in all of 
the stars. 

24 See http://www.sdss3.Org/drl0/irspec/aspcap.php#aspcap 

25 https://www.sdss3.org/drl0/ 

26 MOOG is available at http://www.as.utexas.edu/~chris/moog. 
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Fig. 2.— Measurement of 12 C/ 13 C from 13 C 14 N lines for the 
Li-rich star (circles). The wavelength intervals of S13 are marked. 
The synthetic spectra (colored lines) span a range of 12 C/ 13 C. Top: 
Low 12 C/ 13 C are ruled out in the S13 wavelength interval, but 
continuum level uncertainties make it difficult to place a stringent 
lower limit. A neighboring region (arrow) suggests 12 C/ 13 C~20. 
Bottom: The wavelength interval of S13 favors 12 C/ 13 C~ 30. 


To measure 12 C/ 13 C, we only used the two 13 C 14 N 
lines in S13 because the 13 C 16 0 features are weak and 
contaminated with night sky emission lines. The 13 C 14 N 
lines are also quite weak in all of the stars analyzed, 
which indicates high 12 C/ 13 C. We illustrate the measure¬ 
ment of 12 C/ 13 C from the Li-rich star’s spectrum in Fig¬ 
ure 2. Low 12 C/ 13 C (< 15) are clearly excluded by the 
weakness of the 13 C 14 N. The fits favor 12 C/ 13 C~20-30, 
consistent with normal levels of isotopic 13 C enrichment 
in the stellar envelope following FDU. The constraints 
against high 12 C/ 13 C, i.e., > 30, are obviously weaker 
than those against low 12 C/ 13 C, but our C/N measure¬ 
ment gives us confidence that FDU has occurred. 

In Figure 3, we compare our mixing indicators (C/N 
and 12 C/ 13 C) to red giants in other open clusters. We 
plot the mixing indicators as a function of the clusters’ 
turn off mass (Mto)- The literature values come from 
a series of papers focusing on the abundances of open 
cluster stars (Mikolaitis et al. 2011a,b, 2012). The first 
hfmlpaper compiles even earlier studies from Gilroy (1989); 
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Fig. 3.— Relationship between cluster turn-off mass and C/N 
(top) and 12 C/ 13 C (bottom). The Li-rich star (asterisks) and com¬ 
parison stars (squares) are compared to other open clusters in the 
literature (Mikolaitis et al. 2011a, circles; 2011b, triangles; 2012, 
X). 

Luck (1994); Tautvaisiene et al. (2000); Tautvaisiene et 
al. (2005); Smiljanic et al. (2009) and Mikolaitis et al. 
(2010). Figure 3 demonstrates that the stars in this pa¬ 
per generally fit within the trends defined by the other 
clusters. The exception is C/N for the Li-rich star, which 
is somewhat high for a star of its presumed mass. Oxygen 
is not processed in stars of this mass range, so we expect 
that [O/Fe] to remain near zero. We find that the com¬ 
parison stars have [O/Fe] between —0.03 and 0.11, which 
is consistent with our expectation given the uncertainties 
of ~ 0.1 dex. However, the Li-rich star again stands out. 
It has [0/Fe]= 0.3 ± 0.16 dex, suggestive of an unusual 
enhancement in O. Extra deep mixing would be expected 
to lower the [O/Fe], e.g., as in the case of IRS 7 (Carr et 
al. 2000). 

3.2. Na, Mg, Al, Si, Ti, Fe, and Ni 

The full line lists in Table 5 of S13 were used for Fe I, 
Mg I, and Al I. Subsets of the S13 list were used for Si I 
(15376.8, 15960.1, 16060.0, 16094.8, 16215.7, 16680.8, 
16282.2 A), Ni I (15632.7, 16584.4, 16589.3, 16673.7 A), 
and Ti I (15543.756 A and 15602.842 A). Sodium abun¬ 
dances were measured from two Na I lines (16373.9 and 
16388.9 A) identified in Cunha et al. (2014). We fixed the 



Fig. 4.— Elemental abundances normalized to solar (Asplund et 
al. 2009) of the Li-rich star and three comparison stars. 

C, N, and O abundances found in Section 3.1 for these 
analyses. The abundance results are given in Table 1 
for the main sample, and the uncertainties are again the 
standard deviation of the measurements. 

In Table 2 we compare our results for 
2M19411476+4011008 to the abundances derived 
from optical data by B01. Columns two and three list 
the abundances and standard deviations (a) measured 
here and in B01, respectively. When only one line 
for a given element was measured by B01, we do not 
list a. The last column shows the differences between 
the abundances derived here and those derived by 
B01. Most of our abundances agree within the quoted 
uncertainties. The exceptions are Si and Na. For both 
of these elements, we find lower abundances than B01. 
Compared to Comps 1-3, we generally measure slightly 
larger abundances for 2M19411476+4011008, although 
only Al disagrees with the other cluster stars outside 
the uncertainties. 

In Figure 4, we present our elemental abundance mea¬ 
surements (solar values from Asplund et al. 2009) of the 
Li-rich star and the three comparison stars. Since all of 
these stars have similar evolutionary stages and stellar 
parameters, we would expect them to all have similar 
abundances. (We exclude the star overlapping B01 from 
this plot because it is at a much more evolved evolu¬ 
tionary stage.) The abundances of the Li-rich star are 
generally in good agreement with the Li-normal compar¬ 
ison stars, with the exceptions being O, Na, and Al. The 
enhanced A(O) was already mentioned in the previous 
section. The A(A1) of the Li-ricli star is ^0.14 dex lower 
than the three most similar comparison stars—well out¬ 
side the 0.05 dex uncertainties. A(Na) is also quite low 
in the Li-rich star, though the uncertainties are large. 
However, this large uncertainty means that the two Na 
lines disagree. If we directly compare each line in the 
Li-rich star to the most similar control star, we find that 
both Na lines are significantly weaker in the Li-rich star, 
pointing to a true deficiency in Na. Thus, both A(A1) 
and A(Na) may be low for the Li-rich star. If the reverse 
were true, the anomalous abundances could be explained 
by extra-mixing, which enhances the surface abundance 
of Li, Na, and Al. 
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3.3. Abundance Sensitivities 

To estimate the sensitivity of the abundance measure¬ 
ments to uncertainties in the stellar parameters used to 
select a model atmosphere for each star, we selected 
a fiducial model of T e g = 4700 K, logg = 2.7 dex, 
[M/H]=0.0, and £ = 1.5 km s -1 and tested how the 
output abundances change for perturbations of AT e ff = 
±80 K, A log (7 = ±0.2 dex, A[M/H] = ±0.1 dex, and 
A£ = ±0.2 km s -1 . This was accomplished by adopting 
mean abundances of each element from our stellar sam¬ 
ple, calculating the equivalent width needed to produce 
that abundance with the fiducial model for each line, and 
fixing the equivalent width to see how the output abun¬ 
dance changed with changes in the model atmosphere. 
The changes for all lines of a given element are then av¬ 
eraged and presented in Table 3. Each entry in the table 
gives two numbers. The first number is the change in 
abundance when increasing the stellar parameter by the 
specified amount, and the second number is the change in 
abundance when the stellar parameter is decreased. For 
the sensitivities of C, N, and O, we adopted a simpli¬ 
fied approach by using the line sensitivities of the three 
molecules present in our line list as a proxy for either C, 
N, or O. The molecules in Table 3 have their associated 
elements given in parentheses. 

4. SURFACE GRAVITY 

4.1. Spectroscopic log g 

For elements with low ionization potential that are pre¬ 
dominantly ionized in the photospheres of the stars of in¬ 
terest, the strengths of lines of singy ionized atomic lines 
are sensitive to the log g of the star, while neutral atomic 
lines are not. When lines of both species of an element 
are available, they can be used to measure log g. In the 
APOGEE spectral range, the Ti II line at 15873.834 A is 
the sole known ionized line. Wood et al. (2014) reported 
an experimentally measured log gf = —1.90 ± 0.08 for 
this line, while the APOGEE team adopted an astro- 
physical log gf = —2.06 from fitting spectra of the Sun 
and Arcturus. The latter is adopted here. We measured 
the strength of this Ti II feature and the Ti I lines for our 
five stars while varying the log g of the atmosphere mod¬ 
els from 2.0 to 3.5 dex, as illustrated in Figure 5 for the 
Li-rich star. Our spectroscopic log g is the one for which 
A(Ti) of Ti I and Ti II are equal. These log g measure¬ 
ments are listed in the second section of Table 1, and 
they generally agree well with those predicted from the 
cluster isochrone. The exception is the Li-rich star. We 
find log g = 2.35 dex compared to the isochrone value 
of 2.8 dex; however, this is a good match to the logg 
derived from asteroseismology (Section 4.2). 

To check the impact of having a single ionized line with 
an uncertain oscillator strength, we check for a system¬ 
atic offset in log <7 using the star that overlaps the B01 
study. They measured log <7 with the ionization balance 
of Fe and find log g = 2.6 dex. This agrees with both our 
ionization balance (2.58 dex) and isochrone (2.6 dex) re¬ 
sults. 

4.2. Asteroseismic log g 

An additional gravity measure comes from the astero¬ 
seismic I'm ax in the APOKASC catalog (Pinsonneault et 
al. 2014), using the standard scaling relationship where 
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Fig. 5.— Ionization equilibrium plot for the Li-rich star. The 
shaded regions around Ti I (circles) and Ti II (triangles) show in¬ 
dividual line uncertainty measurements of ±0.05 and ±0.03 dex, 
respectively. The intersection of these relations (x) yields the log g. 
Vertical lines mark the expected log g from isochrones (dotted line) 
and asteroseismology (dash-dotted line). 

g oc v maxvT'ff (Kjcldsen & Bedding 1995). These 
gravities are listed in the bottom section of Table 1. 
In all cases, the asteroseismic gravities agree well with 
our spectroscopic log < 7 . In turn, the spectroscopic log g 
matches well with the isochrone measurement for all but 
the Li-rich star. The Li-rich star has a spectroscopic 
log (7 that is in good agreement with the asteroseismic 
value, but it is 0.4 dex lower than what is expected from 
isochrones. Unfortunately, 2M19411476±4011008 does 
not have a i^ ma x measurement so we do not have a direct 
comparison between the asteroseismic log g and log g de¬ 
rived from ionization balance of Fe lines. 

5. ROTATIONAL VELOCITY 

While measuring abundances, we required a slightly 
larger broadening of the synthetic spectra to fit the lines 
of the Li-rich spectrum compared to the other spectra. 
We acquired the rotational velocities of the stars in this 
study from the pipeline being developed for APOGEE 
(Bizyaev et al., in prep.). Briefly, the spectra are cross- 
correlated with synthetic templates that have been con¬ 
volved with a rotational broadening profile for a range of 
usini. The results of this analysis are listed in Table 1. 
As expected, the Li-rich star has the largest usini, and 
the comparison stars have v sin i near the detection limit 
for APOGEE’s resolution. The clear measurement of 
usini= 8.5 km s _1 points to unusually fast rotation in 
the Li-ricli star. Such a high usinf for RGs is rare. In a 
recent survey for fast rotators among open cluster RGs, 
Carlberg (2014) found only one RG out of ~ 270 cluster 
members that rotate as rapidly as the Li-rich star. 

6 . BINARY COMPANIONS OR OTHER BLENDS 

We performed multiple checks of whether our spec¬ 
troscopic analysis is affected by either a companion to 
the Li-rich star or an unrelated fore/background object. 
First, we fit model spectral energy distributions (SED) 
to photometry spanning optical (Hole et al. 2009), near- 
IR (2MASS, Skrutskie et al. 2006), and mid-IR (WISE, 
Wright et al. 2010) wavelengths, allowing the tempera- 
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TABLE 3 

Abundance sensitivities for the fiducial model 


Species 

AA/A T eff 
±80 K 

AA/A log g 
± 0.2 dex 

AA/A[M/H] 

±0.1 dex 

AA/A$ 

±0.2 km s -1 

CO (C) 

+0.049/—0.055 

+0.065/—0.063 

+0.029/—0.042 

—0.003/+0.002 

CN (N) 

+0.037/—0.021 

+0.053/—0.035 

+0.017/—0.033 

—0.014/+0.015 

OH (O) 

+0.101/—0.097 

—0.017/+0.040 

+0.087/—0.107 

—0.002/+0.002 

Nal 

+0.042/—0.042 

—0.007/+0.027 

—0.002/—0.016 

—0.002/+0.003 

Mg I 

+0.042/—0.038 

-0.041/+0.064 

+0.018/—0.036 

—0.012/+0.013 

All 

+0.070/—0.071 

-0.074/+0.102 

+0.020/—0.044 

—0.033/+0.030 

Si I 

+0.021/—0.015 

-0.044/+0.062 

+0.040/—0.053 

—0.024/+0.023 

Ti I 

+0.111/—0.109 

—0.001/+0.016 

—0.005/—0.012 

—0.045/+0.054 

Ti II 

—0.019/+0.022 

+0.078/—0.092 

+0.035/—0.036 

—0.035/+0.038 

Fe I 

+0.050/—0.045 

—0.032/+0.046 

+0.026/—0.042 

—0.040/+0.043 

Ni I 

+0.019/—0.014 

+0.006/+0.008 

+0.020/—0.033 

—0.017/+0.018 


ture of the models and visual extinction to be free pa¬ 
rameters. The best results are obtained for a bare pho¬ 
tosphere with a hotter temperature (T e 5000 K) than 
that derived from our (J — K s )o calibration and a larger 
line-of-sight extinction. This solution implies that the 
Li-rich star is a more distant star affected by additional 
extinction. However, the analysis favors the largest ex¬ 
tinction we allowed in the fit, Ay = 0.665, which is 
slightly higher than Ay = 0.623, the extinction in the 
Schlegel et al. (1998) map at this position. Additionally, 
the same SED analysis of the most similar comparison 
star also favors a hotter, more reddened solution. There¬ 
fore, we conclude that the temperature difference may 
be systematic and cannot discriminate the Li-rich stars’s 
cluster membership. We also find no evidence for an IR 
excess out to 10 /jrn. 

APOGEE made three unique observations of the Li- 
ricli star, roughly evenly spaced over the course of a 
month. The RV is constant within the 0.1-0.2 km s^ 1 
uncertainties of the three observations, and the Li-rich 
star’s RV matches that of the cluster, as illustrated in 
Figure 1. Admittedly, one month is a short time span to 
look for RV variations, but the literature contains addi¬ 
tional RV measurements for this star spanning at least 
20 years. The optical spectrum taken by AT13 was taken 
between 2010-2012 and also had an RV consistent with 
cluster membership. Hole et al. (2009) define the star as 
a single, cluster member based on four RV measurements 
going back at least to 1992. Taken together, all of these 
RVs suggest that 2M19411367+4003382 is a RV stable 
cluster member. 

We also leveraged the individual visit spectra to search 
for evidence of a secondary spectrum by cross-correlating 
each visit spectra with a library of template spectra. All 
three give best fits to templates of T e g ~ 5000 K. The 
shape of the cross-correlations peaks do not change be¬ 
tween visits, and no secondary peaks are found, ruling 
out bright companions of similar spectral type. A second 
round cross-correlations used the TODCOR algorithm 
(Zucker & Mazeh 1994) to generate a two-dimensional 
cross-correlation using a wide range of possible secondary 
spectral types: no significant secondary peaks were iden¬ 
tified. We conclude that the spectral absorption features 
are truly arising in the red giant’s atmosphere and are 
not significantly contaminated by flux from any binary 
companion or unassociated blend. 

7. INTERPRETATION OF THE LI-RICH STAR 


Our analysis of the Li-rich star has confirmed that it is 
evolved- its low logg and non-solar C/N are consistent 
with the picture of a RG that has undergone FDU and 
should have depleted Li. However, given the uncertain¬ 
ties in the cluster membership, the Li-rich star may be a 
more massive, more luminous background RG that coin¬ 
cidentally shares the same RV, [Fe/H], and CMD locus as 
the cluster. Therefore, to fully explore the implications 
of our measurements we provide two scenarios. The first 
scenario presumes that the Li-rich star is not a member of 
NGC 6819, and we explore how that assumption changes 
our derived abundances, logg, etc. The second scenario 
explores the implications of the star being a cluster mem¬ 
ber, and we demonstrate why we favor this scenario. 

7.1. Case 1: Cluster Non-Membership 
7.1.1. Underestimated Reddening 

If the Li-ricli star is not a member of the cluster, then 
the temperature derived from the J — K s is a lower limit, 
since the reddening may be underestimated. We can 
test how our analysis will change if we adopt the to¬ 
tal line-of-sight reddening observed in the direction of 
2M19411367+4003382. If we adopt Ay = 0.665 (the 
preferred solution of the SED fitting), then E(B — V) = 
0.215 and E(J — K s ) = 0.107. Reapplying our temper¬ 
ature calibration yields 4840 K. The strong temperature 
dependence of Ti I lines increases A(Ti I), while Ti II re¬ 
mains nearly constant. Thus, the ionization equilibrium 
plot (Figure 5) for T e g- = 4840 K would have a similar 
shape, but the horizontal line would shift upwards by 
~ 0.2 dex and intersect the A(Ti II) line at logg = 2.81 
and A(Ti)= 5.14 dex. This means that if the Li-rich star 
is in indeed a background red giant, then its surface grav¬ 
ity is higher than that predicted by the asteroseismology 
and matches that expected from the stellar isochrones. 

Since the Li-rich star has the same metallicity as the 
confirmed cluster stars, we can adopt solar metallicity 
stellar evolution tracks (Bertelli et al. 2008, 2009) to find 
the range of masses and ages of stars near this T e g and 
logg but that have luminosities exceeding that adopted 
in Table 1 under the assumption of cluster member¬ 
ship. We find that stars within T e g = 4840 ± 80 K, 
logg = 2.81 ± 0.2 dex, and logL/Lg > 1.52 span a mass 
range of 1.8-3.0 M 0 and an age range of 0.4-1.6 Gyr. 
The lowest mass stars in this range must be RC stars 
rather than first ascent stars in order to meet all of the se¬ 
lection criteria. First ascent stars are too faint. At higher 
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masses, both first ascent stars and core He burning stars 
satisfy all three constraints. For the entire mass range, 
the stars should have completed first dilution (see, e.g., 
Fig. 1 of Clrarbonnel & Balacliandran 2000) meaning 
the high Li abundance is still unexplained by standard 
stellar evolution. 

Using the sensitivities calculated in Section 3.2, we can 
conclude that increasing the temperature by 140 K has 
the following effects: Fe changes by +0.09, C by +0.09, 
N by +0.06, O by +0.18, Na by +0.07, Mg by +0.07, Al 
by +0.12, Si by +0.04, Ti (from ionization balance) by 
+0.16, and Ni by +0.03. The new C/N is 2.5, which is 
still consistent with the completion of FDU. 

7.1.2. Internal Li Replenishment 

We can use mixing signatures to explore whether in¬ 
ternal replenishment of Li is a viable mechanism for the 
full range of possible masses and ages of the Li-rich star. 
Constraints on these stellar parameters were made in the 
last section for the high extinction case, but we must also 
consider the case where the Li-rich star is affected by the 
same amount of extinction as the cluster stars, i.e., if 
there is negligible extinction between the cluster and the 
Li-rich star. We repeat the analysis with stellar evolution 
tracks for T e g = 4700 ± 80 K, log g = 2.35 ± 0.2 dex, and 
log L/Lq > 1.52. The resulting mass range is M = 1.8- 
4.0 Mq and an age range of 0.2-1 .6 Gyr. In this scenario, 
both first and second ascent evolutionary stages for the 
full mass range are consistent with the selection criteria. 
As a result, 1.8 Mq stars near the cool end of our T e g 
range (4620 K) are luminosity bump stars that also meet 
the gravity and luminosity constraint. 

In either assumed extinction scenario, the Li-ricli star 
is an evolved giant that should have undergone convec¬ 
tive dilution and show depleted Li. The allowed pa¬ 
rameters of the Li-rich star includes two regimes where 
Li regeneration has been suggested to occur—the lumi¬ 
nosity bump, and the narrow mass range of RC stars 
identified by Kumar et al. (2011). Thus, internal Li re¬ 
generation seems to be a viable explanation. However, 
a common feature of internal Li regeneration models is 
that the Li synthesis must happen deep enough in the 
star that whatever mechanism is ferrying the newly syn¬ 
thesized Li into the convection zone must also ferry 13 C 
and thereby lower the surface 12 C/ 13 C. The absence of 
strong 13 C features excludes the deep mixing that should 
accompany Li-regeneration via the D12 models favored 
by AT13, unless we are observing the star near the begin¬ 
ning of the rapid Li enrichment stage (in which the star 
spends less than a few Myr with enriched Li but FDU 
levels of 12 C/ 13 C). 

7.1.3. External Li Replenishment 

External replenishment models run into difficulties 
with the amount of material required to raise the A(Li) 
to the observed value. By rearranging the equation for 
predicting enhanced Li from planet engulfment in Carl¬ 
berg et al. (2012, their Eq. 2) we can write 

l 0 A(Li) 

new _ 10 A ( Li )* 

^ ~ 10 A ( Li )p - 10 A ( Li )new ’ ^ ^ 

where q e is the ratio of the mass of the polluting mate¬ 
rial to the mass in the stellar convective envelope, and 


the A(Li) subscripts of “new,” *, and “p” refer to the 
post-engulfment, the normal stellar, and pollutant abun¬ 
dances, respectively. Adopting A(Li)* = 1.49 dex (the 
largest non-enriched Li), A(Li) new = 2.3 dex (the ob¬ 
served abundance), and A(Li) p = 3.3 dex, (the solar neb¬ 
ular value) we find that q e = 0.09. We used the MESA 
stellar evolution code (Paxton et al. 2013, version 6794) 
to generate limiting mass models of M = 1.7 and 4.0 Mq 
to determine the amount of mass in the outer convection 
zone for stars near 4700 and 4840 K at different phases of 
the RGB evolution. In all cases, that mass was at least 
1.1 Mq. Consequently, the amount of material required 
to explain the enhanced Li in 2M19411367+4003382 is 
0.1 Mq, which is a stellar mass companion. Such a low 
mass companion would have been missed by our binary 
analysis. However, low mass stars are also Li depleted. 
Substellar companions retain their birth supply of Li, 
but they are too low mass unless they are significantly 
enriched in Li over the solar nebular value. 

7.2. Case 2: Cluster Membership 

The inferences on the origin of the excess Li mentioned 
in the previous section are valid in the case of cluster 
membership, but we place further constraints on the un¬ 
usual nature of the Li-rich star if it is truly a cluster 
member. The cluster membership is supported from a 
number of lines of evidence. Its RV is consistent with 
membership, and its optical and infrared magnitudes co¬ 
incide with the cluster’s CMD locus. This latter piece 
of evidence is further supported by the fact that the 
temperature we derived from the 2MASS colors is the 
same as that derived by the ASPCAP pipeline (see Sec¬ 
tion 2). The ASPCAP pipeline fits model spectra to 
the continuum-normalized observed spectrum and is, im¬ 
portantly, independent of any assumption of reddening. 
Therefore, the ASPCAP temperature verifies that the 
reddening of the Li-rich star is the same as the rest of 
the cluster and thus shares the temperature and magni¬ 
tude of the cluster’s isochrone. 

7.2.1. A Deeper Look at Proper Motions 

One of the major pieces of evidence against member¬ 
ship for the Li-rich star is that Platais et al. (2013) find 
proper motions inconsistent with membership. However, 
the authors of that paper noted that finding low proper 
motion membership probabilities for stars with high RV 
membership probabilities is likely due to source confu¬ 
sion (see their Figure 6 ). In fact, an earlier proper motion 
study (Sanders 1972) gave the Li-rich star (star 90 in that 
paper) a 90% membership probability. An inspection of 
a Digitized Sky Surveys (DSS) optical blue image and 
the 2MASS H -band image in Figure 6 reveals an elonga¬ 
tion in the optical DSS image that is not present in the 
M-band image. This is suggestive of a faint, blue source. 
If this source is visible in the plates used by Platais et al. 
(2013), or worse is only visible in some of the plates, it 
could shift the apparent photocenter of the Li-rich star 
and skew the proper motion measurement. Star 90 on 
the Sanders (1972) plate is very round, suggesting that 
the nearby blue source is too faint to adversely affect 
the earlier proper motion measurement. Indeed, the V 
limit of the Sanders (1972) plate is quoted as ~ 14.5 mag, 
whereas the magnitude limits of the plates used in Platais 
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Fig. 6.— Image of the Li-rich star in 2MASS H -band image (solid 
lines) and the DSS optical blue image (dashed lines). The optical 
image of the Li-rich star is slightly elongated on one side compared 
to the infrared image. 


et al. (2013) are V = 17-22 mag. Because this visual 
companion is faint and blue, it will not contribute sig¬ 
nificantly to the H-band spectrum of the Li-rich star, 
consistent with our analyses in Section 6 . 

7.2.2. Problems with Replenishment Mechanisms 

Cluster membership adds more problems to both the 
internal and external replenishment models. As a clus¬ 
ter member, the evolutionary stage of the Li-rich star 
is clearly below the luminosity bump, which led A13 to 
favor the recent D12 models. However, there is a dis¬ 
crepancy between the D12 models and the Li-rich star’s 
CMD position. D12’s model predicts an evolution that 
traces an extended luminosity bump-like zigzag on the 
CMD. Unlike the luminosity bump, the extended zigzag 
in the D12 models parallels the first ascent RGB at lower 
luminosity. For a given T e g , the A log L is « — 0.2, corre¬ 
sponding to a magnitude difference of +0.5. The Li-rich 
star does not sit 0.5 magnitudes below the RGB in either 
optical or infrared CMDs. For external replenishment 
models, the other cluster RGs provide a revised estimate 
for pre-engulfment A(Li), which is lower than that esti¬ 
mated in Section 7.1.3. Additionally, with a well known 
evolutionary stage, we can more accurately pinpoint a 
MESA model to estimate the size of the convection zone. 
Using the 1.7 Mq MESA run, we inspected the internal 
profiles for four evolutionary stages on the first ascent 
RGB near T = 4700 K. In all cases, there is ~ 1.4 Mq in 
the convective envelope, requiring 0.15 Mq of accreted 
material to explain the observed A(Li). 

7.2.3. An Unusually Low Mass for the Li-rich Star 

Another implication of the cluster membership is that 
it allows us to put stronger credence on the asteroseis- 
mology result, since our spectroscopic log <7 agrees with 
the asteroseismic log < 7 . The scaling relations for ^ max 
and A v can both be cast to depend only on T e g, M, 
and either R or L. Since T e s is accurately known and 
has the weakest influence on zz max and A u, the remaining 
quantities can be derived without any other assumptions. 
These results are given in the bottom section of Table 1. 


There is good agreement between the luminosity derived 
through asteroseismology and that derived assuming the 
star is in the cluster. This analysis places the Li-rich star 
at roughly the distance of the cluster but with a mass of 
only 0.7 Mq, which explains the “low” log < 7 . The main 
sequence lifetime of such a low mass star exceeds the age 
of the universe. If this mass is correct, the most plausible 
explanation is that the star lost a significant fraction of 
its mass, and we suggest that this mass loss is related to 
the source of Li. Furthermore, this low mass does not 
nullify the CMD position as evidence of cluster mem¬ 
bership. The luminosity of stars with H-burning shells 
is dictated by the mass of the stellar core, leading to a 
core mass-luminosity relationship. Boothroyd & Sack- 
mann (1988) demonstrated that this relationship holds 
for intermediate mass stars and that the total mass has a 
negligible effect on this relationship. Thus, if the Li-rich 
star began its life with the same mass as the other RGs 
and only lost mass recently, it should have the same core 
mass as the other RGs. 

7.2.4. External Enrichment at Low Mass 

The case for enrichment by an external source of A(Li) 
is much more feasible if the star is really 0.7 Mq. We 
can again adopt A(Li) <0.7 dex for the pre-engulfment 
abundance and use the MESA models of the 1.7 Mq 
star. If an engulfed companion strips 1 Mq of mate¬ 
rial from the star before it disperses into the remain¬ 
ing stellar envelope, then the remaining envelope is only 
0.4 Mq. Since the star is not entirely stripped of the 
original convective envelope, the mass loss not does ex¬ 
tend deep enough to expose material with low 12 C/ 13 C. 
Only a 45Mj up object is required to enrich A(Li) up to 
the observed value. This mass is safely in the substel- 
lar companion regime of objects, which retain their birth 
A(Li). 

7.2.5. Updated. Abundances at Low log g 

If the 2M19411367+4003382’s log g is truly so low, then 
we must re-derive the abundances. Because using the 
abundance sensitivities in Table 3 requires extrapolat¬ 
ing by 0.2 dex beyond the log g range we tested, we in¬ 
stead calculate all of the new abundances with the lower 
logg atmosphere. We calculate the equivalent widths 
needed to create the abundances we measured for the 
Li-rich star with the log g = 2.80 dex model and then 
find how the abundances changes for these same equiv¬ 
alent widths with a log <7 = 2.35 dex model. For the C, 
N, and O abundances from the molecular lines, we iter¬ 
ate the process until the output abundances match the 
input abundances. The lower log g has the following ef¬ 
fects: Fe changes by +0.09, C by —0.15, N by +0.08, O 
by +0.02, Na by +0.03, Mg by'+0.16, Al by +0.17, Si 
by +0.11, Ti by +0.05, and Ni by —0.005 dex. Figure 7 
shows these new abundances compared to that of the 
comparison stars. The updated A(A1) is now in excel¬ 
lent agreement with the other cluster stars. The rather 
large increases in the abundances of Mg, Si, and Fe are 
still consistent with the abundances of the other cluster 
stars within the uncertainties. Furthermore, the C/N 
ratio is now 1.4, which is also in much better agreement 
with other RGs in both NGC 6819 and in other clusters 
(Figure 3). However, oxygen is still unusually high, with 
[0/Fe]~ 0.23 ± 0.15 dex. 




10 


Carlberg et al. 



Fig. 7.— Elemental abundances normalized to solar (Asplund 
et al. 2009). The Li-rich star’s abundances are computed for a 
T e fj=4700 K and logg=2.35 dex atmosphere model. The abun¬ 
dances of the control stars are the same as in Figure 4. 

7.2.6. Comparison to Other Known Li-rich Cluster Giants 

Finally, if 2M19411367+4003382 is truly a member of 
NGC 6819, then we can compare this star to the few 
other known Li-rich red giants in other open clusters. 
These include Trumpler 5 (TV 5) #3416 (Monaco et al. 
2014) and Berkeley 21 (Be 21) T33 (Hill & Pasquini 
1999). Both of these open clusters have significantly sub 
solar metallicities, near —0.5 dex. In both cases, the Li- 
rich stars are super Li-rich, A(Li)> 3.3 dex. Tr 5 3416 is 
also known to be a slow rotator and have 12 C/ 13 C=14. 
The Li-rich star in Tr 5 sits on the RC, while the 
star in Be 21 sits a magnitude above the RC. There¬ 
fore, 2M19411367+4003382 does not share many com¬ 
mon properties with these other known Li-rich stars in 
open clusters. It is solar metallicity, fainter than its host 
cluster’s RC, not super Li-rich, and has faster rotation 
and higher 12 C/ 13 C than the Tr 5 Li-rich star. Addition¬ 
ally, 2M19411367+4003382’s initial mass was presum¬ 
ably 1.7 Mq, and Monaco et al. (2014) estimates that 
the masses of the Tr 5 RC stars are 1.0-1.4 Mq, while 
Bragaglia & Tosi (2006) gives 1.4 Mq as the mass of the 
turn-off stars in Be 21. 

8 . CONCLUSIONS 

We have performed a spectroscopic analysis of the 
Li-rich star 2M19411367+4003382, and we have pro¬ 
vided crucial measurements of log g, C/N, and 12 C/ 13 C 
that confirm that the star is an evolved red giant with 
evidence of having completed FDU (in its C/N and 
12 C/ 13 C) but showing no evidence of low 12 C/ 13 C that 
should accompany the previously favored model of ex¬ 
plaining the star’s enriched Li. We have also measured 
abundances of eight additional elements and found a 
chemical similarity between the Li-rich star and the other 
cluster members in elements that are not affected by nu¬ 
clear processing, including both iron group (Fe, Ni) and 
a elements (Mg, Si, Ti). 

The membership of 2M19411367+4003382 has recently 
been questioned by both asteroseismic and recent proper 
motion measurements. We have argued that an older 
proper motion measurement may be more reliable, and 


this older measurement is compatible with member¬ 
ship. The asteroseismic non-membership evidence ac¬ 
tually contains a surprising result. The asteroseismic 
i/max yields a log <7 lower than that expected for mem¬ 
bers, which we confirm spectroscopically. Recomputing 
the abundances at lower log g brings the abundances of 
the Li-rich star in better agreement with the other red 
giants in NGC 6819. Furthermore, v max and Av also 
gives the star’s luminosity and mass and places the Li- 
rich star at roughly the distance to the cluster, but with 
a mass of only 0.7 Mq. Such a low mass star could not 
have evolved to the RG stage, implying that the Li-rich 
star has lost significant mass. The similar luminosity to 
other NGC 6819 RGs implies that the Li-rich star has a 
similar core mass to the other RGs, which is also con¬ 
sistent with the picture of extensive mass loss. For the 
case of external enrichment, this low stellar mass is the 
only case where a putative engulfed companion can have 
low enough mass to have retained its birth Li but still 
be massive enough to contribute enough Li to explain 
the observed stellar A(Li). A more detailed analysis of 
the stellar oscillation spectrum may provide new insights 
on whether the scaling relationships do not apply to this 
star or whether the internal structure of the Li-rich star 
can shed light on its anomalous Li abundance and mass. 
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